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Abstract 


This  is  the  final  technical  report  on  ONR  Grant  N00014-89-J-1794,  entitled  "Application  of  Fractal 
Concepts  to  the  Analysis  of  Atmospheric  Processes",  for  the  period  from  1  April  1989  through  30 
September  1991.  The  goal  of  the  research  was  to  improve  understanding  of  atmospheric  processes 
through  the  application  of  analytical  approaches  derived  from  methods  that  have  been  associated 
with  fractal  and  wavelet  methods.  This  was  the  first  phase  of  ongoing  work  that  is  now  supported 
by  a  new  grant  for  the  period  to  30  September  1994. 

There  were  two  major  thrusts  of  this  work  to  date.  First,  we  have  applied  multiresolution  feature 
analysis  to  an  atmospheric  data  set  In  the  process  the  analysis  was  extended  to  three-dimensional 
vector  fields.  To  m^e  the  three-dimensional  vector  features  "imaginable,"  we  have  adopted  the  a 
variant  of  the  Lorenz  concept  of  empirical  orthogonal  functions  or  EOF.  These  allow 
representation  of  empirically-defined  patterns  of  small-scale  atmospheric  motions.  Second,  we 
developed  a  strategy  for  the  generalized  wavelet  analysis  of  homogeneous  turbulence.  The  strategy 
calls  for  decomposition  of  a  collection  of  direct  numerical  simulation  data  sets  for  homogeneous 
turbulence  in  a  stratified  Boussinesq  fluid  subject  to  an  imposed  shear,  in  both  two  and  tiiree 
dimensions. 
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1 .  Introduction 


This  final  technical  report  for  Grant  N00014-89-J- 1794,  entitled  "Application  of 
Fractal  Concepts  to  the  Analysis  of  Atmospheric  Processes",  covers  the  period  from  1 
April  1989  through  30  September  1991.  Tlie  project  began  as  one  focused  on  fractal 
methods,  but  it  was  expanded  during  the  course  of  the  work.  The  work  described  herein 
represents  the  initial  phases  of  the  project  and  its  status  through  the  end  of  the  above  grant; 
a  new  grant  was  awarded  to  continue  the  work  until  September  1994. 

The  goal  of  the  research  is  to  improve  our  understanding  of  atmospheric  processes  through 
the  application  of  analytical  approaches  derived  from  methods  that  have  b^n  associated 
with  fractal  and  wavelet  analyses.  We  seek,  first,  a  clear  understanding  of  the  spatial  and 
temporal  characteristics  of  atmospheric  motions  that  contain  many  scales  and  are 
intermittent;  second,  we  seek  to  determine  the  causes  of  key  features  of  these  motions. 

Our  near  term  objectives  included: 

Thrust  1;  Multiresolution  Feature  Analysis:  To  identify  and  obtain  existing  data  that  are 
suitable  for  studying  spatial  distributions  of  atmospheric  properties  [including  aerosols 
and  plmne  materials],  to  develop  analysis  methods  and  the  computer  algorithms  and  codes 
necessary  to  apply  them,  to  apply  the  analysis  methods  to  atmospheric  observations,  to 
interpret  the  results  of  the  an^yses  to  identfy  and  characterize  important  physical 
processes,  and  to  communicate  the  results  of  the  analyses  and  interpretations. 

Thrust  2:  Wavelet  Analyses:  To  identify  and/or  create  appropriate  wavelets  for  application 
in  the  analysis  of  field  data  and  direct  numerical  simulations  of  atmospheric-type  flows.  A 
key  question  to  be  answered  [which  has  not  been  rigorously  answered  yet]  is  "Given  a 
signal  or  data  sample,  what  analyzing  wavelet  ought  to  be  used  to  examine  it?"  We  seek  a 
criterion  or  criteria  of  optimality  driven  by  the  objectives  of  the  analysis,  i.e.,  what  do  we 
seek  to  learn? 

Appendixes  1  and  2  contain,  respectively,  a  publication  and  a  list  of  presentations  on  the 
work  supported  by  the  grant.  Appendix  3  contains  a  few  relevant  references. 

2.  Multiresolution  Feature  Analysis 

The  multiresolution  feature  analysis  [MFA]  methodology  developed  by  Jones,  et  al.  (1988) 
was  programmed  and  applied  to  lidar  and  transmissometer  smoke  plume  imagery.  Other 
fracM  dimension  calculation  methods  are  also  being  programmed  and  used  on  these  data. 
The  MFA  method  successively  degrades  image  resolution  by  half,  applies  a  "feature¬ 
detecting  filter"  to  the  degraded  images,  and  examines  the  changes  in  the  number  of 
features  whose  "intensity"  exceeds  different  thresholds.  The  statistical  relationships  among 
features  at  different  scales  define  scaling  parameters  and  fractal  dimensions. 

We  extended  the  MFA  concepts  so  that  they  can  be  applied  to  atmospheric  vector  fields.  In 
essence,  the  usual  scalar  filtering  techniques  have  been  redefined  to  replace  the  arithmetic 
products  of  numbers  with  inner  (dot)  products  between  vectors.  We  had  intended  to  define 
features  that  have  some  special  physical  significance,  but  had  difficulty  speci^ng  more 
than  a  very  few,  so  we  sought  an  objective  method  for  defining  "physically  significant 
features." 

I.x)renz  (1956)  used  a  variant  of  principal  component  analysis  to  define  empirical 
orthogonal  functions  (EOF)  that  represented  different  patterns  of  variability  in  atmospheric 
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pressure  fields  over  the  U.  S.  Ludwig  and  Byrd  (1980)  adapted  this  technique  to  identify 
vector  patterns  that  accounted  for  the  most  variance  in  a  data  set  The  use  of  these 
techniques  for  defining  those  small  scale  patterns  of  variability  that  explain  most  of  the 
variance  seemed  a  logical  approach  to  the  definition  of  features  to  be  used  with  the  MRF 
analysis.  The  program  necessary  to  calculate  the  most  important  EOFs  for  small  sections 
of  three  dimensional  vector  fields  has  been  written  and  applied  to  Schneider's  wind  data 
[see  Appendix  1].  A  corresponding  methodology  for  sc^ar  fields  has  been  programmed 
and  is  being  tested.  Available  lidar  and  transmissometer  smoke  plume  observation  data 
appear  to  have  artifacts  that  lead  to  anomalous  results. 

The  first  task  of  this  study  was  to  identify  and  acquire  suitable  data  for  analysis.  This  task 
has  largely  been  completed.  However,  this  is  the  kind  of  task  that  never  completely  ends 
and  we  will  continue  trying  to  identify  suitable  data  for  analysis.  In  fact,  we  recendy 
identified  a  large  eddy  simulation  (Costigan,  1992)  of  the  same  situation  described  by 
Schneider's  wind  data.  These  results  are  being  obtained.  The  data  acquired  include  large 
eddy  and  direct  numerical  simulation  results,  dual  Doppler  radar  observations  of 
atmospheric  motions,  lidar  cross  sections  and  transmissometer  imagery  through  smoke 
plumes,  and  a  digitized  weather  satellite  image  of  marine  clouds. 

Much  of  the  second  task  which  was  to  identify  analysis  approaches  has  been  done,  but  we 
continue  to  identify  potentially  useful  approaches  that  deserve  to  be  pursued.  Computer 
programs  have  been  written  for  calculating  the  fractal  dimension  of  two  dimension^  scalar 
fields  by  box  counting,  Fourier  filtering  and  multiresolution  feature  analysis  methods.  The 
programs  have  also  been  written  for  identifying  physically  significant  features  for  the  MFA 
meAodology.  Some  of  these  programs  have  been  applied  to  various  data  sets  in 
accordance  with  the  third  task,  but  little  of  the  interpretation  task  has  been  done. 

The  primary  product  of  this  thrust  area  to  date  is  the  paper  displayed  in  Appendix  1. 

3 .  Wavelet  Analyses 

To  analyze  scale  interactions  in  a  fluid  dynamical  field  one  would  like  to  have  the  ability  to 
consider  space  and  time  and  their  respective  scales  simultaneously.  This  suggests  using  a 
multiresolution  analysis  based  on  the  wavelet  transform,  an  operation  originally  defined  in 
a  paper  by  Grossman  and  Morlet  (1984).  Our  initial  approach  was  direct,  namely,  to  study 
the  known  and  successful  wavelets,  e.g.,  Morlet,  Mexican-hat,  etc.,  and  to  ascertain  their 
properties  in  a  systematic  way,  with  particular  focus  on  what  the  wavelet-transformed  data 
signal  tells  us.  In  a  majority  of  the  published  reports,  wavelets  have  been  applied  in  an  ad 
hoc  fashion,  with  little  commentary  regarding  the  preference  of  one  wavelet  over  another. 
Our  experiments  with  a  number  of  wavelets  indicated  that  the  choice  of  wavelet  is  fairly 
inconsequential.  Accordingly,  for  our  ongoing  work,  we  have  opted  to  use  a  class  of 
scale-based  filters  which  are  very  nearly  true  wavelets,  namely  Fourier  bandpass  (FBP) 
filters.  In  physical  space  these  fUters  are  the  difference  of  two  appropriately  scaled  sine 
functions  of  different  amplitudes  and  widths. 

Our  strategy  calls  for  the  use  of  wavelets  to  offer  new  insights  in  regard  to  the  dynamics  of 
stratified  flows.  The  data  for  this  analysis  will  be  a  collection  of  direct  numericad 
simulation  data  sets.  We  will  consider  homogeneous  turbulence  in  a  stratified  Boussinesq 
fluid  subject  to  an  imposed  shear,  in  both  two  and  three  dimensions.  This  is  a 
generalization  which  includes  isotropic  decay,  decay  of  turbulence  in  a  stratified  fluid,  and 
shear  of  a  uniform  density  fluid  as  special  cases.  This  scenario  covers  most  cases  of 
practical  interest,  especially  in  the  context  of  smaller  scale  geophysical  flows.  We  envision 
the  application  of  the  wavelet  methodology  to  the  following:  (1)  classification  of  flow 


structure,  (2)  mixing  efficiency  measures,  and  (3)  internal  wave  /  turbulence 
discrimination. 


Appendix  1.  "Subkilometer  Spatial  Variability  of  Winds  in  a  Sheared, 
Convective  Atmospheric  Boundary  Layer"  by  Ludwig,  Schneider  and 
Street 
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SUBKILOMETER  SPATIAL  VARIABILTY  OF  WINDS  IN  A  SHEARED, 
7.3  CONVECTIVE  ATMOSPHERIC  BOUNDARY  LAYER 

by  F.  L.  Ludwig*  ,  J.  M.  Schneider**  and  R.  L.  Street* 


1.  INTRODUCTION 

This  work  evolved  from  an  effort  to  use  fractal  concepts 
to  interpret  atmospheric  behavior  into  one  where  some  of  the 
techniques  used  to  calculate  fractal  dimension  have  been  used 
to  study  the  nature  of  the  air  motions  themselves,  providing 
new  insights  into  Richardson’s  (1922)  little  verse: 

"Big  whorls  have  little  whorls. 

Which  feed  on  their  velocity; 

And  little  whorls  have  lesser  whorls, 

And  so  on  to  viscosity  ..." 

We  have  attempted  to  define  Richardson’s  whorls.  This 
would  not  have  been  possible  without  the  dual  Doppler  radar 
wind  data  from  the  1984  Phoenix  II  project.^ 

2.  DUAL  DOPPLER  RADAR  WIND  DATA 

The  Phoenix  n  data  have  unusual  three-dimensional 
resolution  on  a  grid  with  200  m  spacing  between  points  (in  all 
three  direcuons)  over  a  9x9  km  area  to  a  depth  of  2  to  3  km. 
Temporal  resolution  was  about  two  minutes.  Three-dimen¬ 
sional  wind  vectors  (vertical  components  were  inferred  from 
continuity  assumptions)  are  available.  Some  data  are  missing 
because  the  radar  volumes  did  not  overlap  at  the  comers,  or 
there  was  often  insufficient  chaff  at  the  upper  levels,  but  the 
spatial  resolution  is  still  sufficient  that  we  were  able  to  de¬ 
velop  the  analytical  approaches  discussed  later. 

The  observation  area  (about  20  km  east  of  the  foothills 
of  the  Rocky  Mountain  Front  Range  and  about  50  km  east  of 
a  reasonably  uniform  ridge  of  the  continental  divide)  is  rela¬ 
tively  flat,  with  undulations  of  about  ±50  ra  around  an  aver¬ 
age  elevation  of  16CX)  m  ASL.  The  relative  proximity  to  the 
mountains  appears  to  have  controlled  the  flow.  Weak  east¬ 
erly  surface  winds  (-3  m  s“i)  and  the  strong  isolation  sug¬ 
gested  a  typical  convective  planetary  boundary  layer  (PBL), 
well-mix^  and  capped  by  a  strong  inversion.  Aircraft  and 
radar  data  indicated  otherwise.  Each  of  the  days  (17  and  22 
June  1984)  chosen  for  analysis  had  a  multi-layered  structure. 
The  locally  generated  convective  layer  was  interacting  with  an 
over-running,  slightly  stable  layer  of  moderate  turbulence  that 
had  probably  formed  as  a  convective  PBL  over  the  moun¬ 
tains.  The  westerlies  that  advcctcd  this  over-running  layer, 
combined  with  the  underlying  easterlies  to  maintain  a  strong 
shear  throughout  the  depth  of  the  PBL.  This  shear  persisted 
in  spite  of  the  afternoon  convection  which  suggests  a  ther¬ 
mally  driven  easterly  flow  in  the  lower  layers,  induced  by  the 
warm,  cast-facing  mountain  slopes.  The  strong  persistent 
shear  induced  important  vertical  motions.  The  turbulent  ki- 
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netic  energy  budget  for  these  flows  (Schneider,  1991)  shows 
shear  generation  that  was  frequently  as  great,  or  greater  than 
buoyant  production  above  the  lowest  kiiomcicr.  There  were 
more  than  the  usual  thennais  ascending  until  their  buoyancy 
is  lost  to  entrainment,  or  they  are  stopped  by  an  inversion 
layer  aloft.  Buoyant  forcing  and  shear-generated  eddies  de- 
tennined  vertical  motion. 

The  most  unstable  of  the  periods  analyzed  (22  June 
1984, 12:40-13:(X)  Mountain  Standard  Time)  had  a  rich  mix 
of  vorticity  structures  (Schneider,  1991)  that  evolved  and  in¬ 
teracted  throughout  the  period.  The  largest  and  longest-lived 
of  the  vortex  filaments  were  near  the  level  where  the  mean 
winds  changed  ftom  easterly  to  westerly.  They  were  kilome¬ 
ters  long  and  persisted  for  more  than  ten  minutes,  with  an 
orientation  that  was  roughly  horizontal  and  normal  to  the 
shear,  i.c.  approximately  north-south.  Shorter  (in  both 
length  and  duration)  vorticity  features  parallel  to  the  shear,  at 
altitudes  above  and  below  the  major  features  also  appeared. 

There  were  more  stable  cases  too,  when  the  convection 
was  much  less  vigorous.  These  boundary  layers  were  decay¬ 
ing  because  of  subsidence,  or  decreased  surface  heat  flux.  In 
one  (17  June  1984, 15:00-15:21  MST),  there  was  very  little 
vertical  motion  or  turbulent  kinetic  energy  production.  That 
case  provides  an  interesting  contrast  with  the  unstable  case 
just  described.  We  focus  on  these  two,  very  different  cases. 

Schneider  (1991)  has  described  the  processing  of  the 
original  data  in  detail.  Each  volume  scan  provided  informa¬ 
tion  from  about  lO^  volumes,  each  about  2.7x10^  (-140 

m  on  a  side).  The  processing  to  obtain  winds  in  Cartesian 
coordinates  inherently  smoothed  the  measurements. 
Cressman  weighted  averaging  was  used  to  get  the  two  radial 
wind  components  (one  from  each  radar)  on  a  Cartesian  grid; 
it  provides  a  nearly  uniform  fUicring  that  is  consistent  with 
the  volume  averaging  introduced  by  the  radar  observations. 
Additional  Leisc  (1981)  filtering  was  used  for  numerical 
smoothness.  Radial  winds  were  linearly  interpolated  between 
successive  observations  to  common  times.  A  standard  pack¬ 
age  of  NCAR  programs  (CEDRIC;  Mohr  et  al.,  1986)  was 
used  to  calculate  u,  v,  and  w  components  from  the  radar  ob¬ 
servations  and  the  continuity  equation. 

3.  MULTIRESOLUTION  FEATURE  ANALYSIS 

3.1  Scalar  Fields 

Jones  et  al.  (1991)  define  multiresolution  feature  analy¬ 
sis  (MFA)  to  have  the  following  steps: 

1 .  Repeatedly  smooth  and  resample  a  scalar  field,  giving  a 

new  field  at  each  step  with  hif  the  resolution. 

2.  Apply  a  feature  detector  to  the  fields  from  step  1. 

3.  Identify  peaks  in  the  fields  produced  in  step  2, 

4.  Count  the  peaks  from  step  3  that  exceed  each  of  a  range 

of  threshold  values  for  each  resolution. 

5.  Rescale  and  plot  the  results  —  (Si)2(ny)i  versus  y/(Si)*'. 

where  Si  is  the  i^h  spatial  scale;  (ny)i  is  the  number  of  peaks 
exceeding  the  threshold  y  at  the  resolution.  A  scaling  field 
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has  a  value  of  k  that  pnxiuces  a  set  of  invariant  relationships; 
k  is  related  to  the  fra^  properties  of  the  distribution  of  the 
feamres.  Our  interest  here  is  in  the  smoothing,  and  the  defi¬ 
nition  and  application  of  feature  detectors,  the  first  two  steps 
in  the  above  list. 

The  smoothed  value  at  a  point  is  obtained  tirom  a 
weighted  average  of  values  at  surrounding  points.  The 
smoothed  field  will  have  cxie-quarter  as  many  points,  each 
representing  four  times  the  area  of  a  cell  in  the  original, 
higher  resolution  field.  The  sinoothing  cliiiiiiiatcs  sriialler 
scale  features.  Feature  detection  as  described  by  Jones  et  al. 
(1991)  works  in  much  the  same  way  as  the  smoothing,  ex¬ 
cept  that  the  value  at  a  point  is  given  by  the  sum  of  the  prod¬ 
ucts  of  the  surrounding  values  with  a  template  of  numbers 
representing  the  feature  to  be  detected.  For  example  a  detec¬ 
tor  for  a  vertically  oriented  "edge"  might  consist  of  a  3x3  ar¬ 
ray  of  numbers,  with  the  first  column  made  up  of  positive 
values,  the  middle  of  zeroes  and  the  third  of  negative  values. 
When  this  ten^late  is  placed  over  a  region  of  strong  gradient 
from  left  to  right,  it  will  produce  a  large  positive  value.  Other 
configurations  can  be  devised  to  respond  to  other  specific 
patterns. 

3.2  Extension  to  Three-Dimensional  Vector  Fields 

Smoothing  is  much  the  same  in  cither  two  or  three  di¬ 
mensions,  differing  for  the  latter  only  in  that  the  filter  will  be 
three-dimensional.  The  definition  of  vector  "features”  is 
somewhat  subtler.  Either  one  can  apply  the  technique  to 
some  scalar  property  of  the  vector  field  (e.g.  divergence  or 
vertical  vorticity  component ),  or  define  vector  features  and 
apply  the  multiresolution  methodology  directly,  i.e.  the  def¬ 
inition  of  features  as  three-dimensional  arrays  of  vectors  over 
small  subregions.  The  latter  approach  was  taken.  The  vector 
feature  detectors  are  defined  so  that  their  output  is  the  sum  of 
the  scalar  products  of  their  component  vectors  with  the  corre¬ 
sponding  vectors  in  the  field  being  analyzed.  The  result  is  a 
three  dimensional  array  of  scalar  "feature  intensities."  One  of 
the  important  differences  between  this  approach  and  conver¬ 
sion  to  scalars  before  the  feature  analysis,  is  that  the  results 
can  be  interpreted  as  specific  kinds  of  "whorls." 

3.3  Empirical  Definition  of  Vector  Features 

It  has  been  difficult  to  imagine  three-dimensional  vector 
features,  so  we  devised  a  method  by  which  the  data  would 
determine  important  patterns  of  sm^  scale  variability.  Our 
approach  goes  by  many  names  (e.g.  principal  component, 
eigenvector,  or  faaor  analysis).  It  represents  the  data  in 
terms  of  a  set  of  linearly  independent  basis  functions,  which 
Lorenz  (1956)  called  empirical  orthogonal  functions  -  EOF. 
This  abbreviation  is  generally  understood,  and  can  be  reinter¬ 
preted  to  mean  empirical  orthogonal  features^  a  more  apt  de¬ 
scription  for  this  application. 

EOF  features  are:  (1)  based  on  the  characteristics  of  the 
data  themselves,  (2)  linearly  independent  of  one  another,  and 
(3)  charaacrized  by  a  measure  of  their  relative  importance. 
Past  EOF  applications  have  focussed  on  scalar  features  or 
sets  of  basis  functions  for  numerical  calculations.  Lorenz’ 
(1956)  patterns  of  atmospheric  pressure  at  64  stations  in  the 
United  States  was  a  scalar  application.  Lumley  (1967,  1980) 
suggested  similar  analyses  to  extract  coherent  structures  from 
turbulent  flows.  He  used  complex  number  representation  of 


vectors  over  the  entire  domain.  Others  (Hardy,  1977; 
Ludwig  and  Byrd,  1980)  used  vector  EOFs  to  classify  inputs 
for  wind  models.  Sirovich  (1988)  analyzed  numerical  simu¬ 
lations  with  EOFs. 


Lorenz  (1956)  provides  a  complete  discussion  of  the  use 
of  EOFs  for  representing  maximum  variance  (i.e.,  informa¬ 
tion)  in  a  data  set  with  fewer  descriptOTS.  Parameter  values 
representing  the  state  of  some  system  can  be  represented  by: 
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The  left  hand  side  (LHS)  is  the  set  of  observed  parameters 
minus  their  individual  mean  values  over  all  cases.  The  q's 
are  specific  to  each  subregion,  and  can  be  used  to  replace  the 
original  descriptors.  If  all  the  q's  are  used,  the  observations 
are  represented  exactly.  If  the  sequence  is  truncated,  the  rep¬ 
resentation  is  approximate.  The  y  vectors  are  the  EOFs,  ar¬ 
ranged  in  decrea^g  order  of  explained  variance.  They  can 
be  interpreted  as  patterns  of  variability  about  the  mean.  The 
coefficients  qic  are  the  inner  (scalar)  products  of  the  observa¬ 
tion  column  vector  (LHS  of  Eq.  1)  with  the  corresponding 
unit  EOF  vector  (the  y's  on  the  RHS).  The  y's  arc  the 

eigenvectors  of  the  covariance  matrix  (P-P)  (P-P)^  formed 
by  multiplying  the  data  (relative  to  their  means)  matrix  by  its 
transpose.  The  eigenvector  calculation  is  straightforward 
with  available  FORTRAN  subroutines  such  as  are  provided 
in  Numerical  Recipes  (Press  ci  al.,  1987). 

The  approach  here  differs  from  others  in  that  it  uses 
small  3x3x3  subsections  of  the  observed  wind  field.  The 
middle  vector  is  subtracted  from  each  of  the  surrounding  26 
vectors  to  obtain  a  difference  vector  (aV)  at  each  point  in  the 
subsection: 
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where  ujjic,  Vijk  and  wyi^  arc  vector  components  at  point  ijk  in 
the  m**'  subsection;  222  is  the  center  point.  The  3-dimcnsional 
array  of  AVs  shows  the  local  pattern  of  motion.  If  we  have 
this  array  and  the  vector  V72?  at  the  center  point,  the  sur¬ 
rounding  26  vectors  can  be  reconstructed.  The  deviations 
AV'  about  the  averages  (for  the  same  relative  points  in  all  the 
N  subsections)  are  described  by  [(3x3x3)  -1]  grid  points  x3 
vector  components,  i.e.  78  components  for  each  column 
vector  describe  the  state  of  a  small  subcube.  The  matrix  of 
these  state  vectors  is  multiplied  by  its  transpose  to  give  the 
covariance  matrix  from  which  the  EOFs  arc  found.  Those 
EOFs  that  explain  the  most  variance  in  the  individual  patterns 
can  be  used  as  the  features  in  the  modified  multiresolution 
feature  analysis  methodology.  They  can  also  be  considered 
in  some  sense  to  represent  Richardson’s  whorls. 


Data  from  several  times  during  a  sequence  (generally 
spanning  --20  minutes)  were  combined  so  that  more  than 
150,(XX}  3x3x3  subvolumes  were  typically  available  for  anal¬ 
ysis.  The  data  matrix  was  formed  horn  a  randomly  selected 
sample  of  -10  percent  of  the  total  to  reduce  confutation  and 
memory  requirements.  We  repeated  the  process  for  ten  dif¬ 
ferent  samples  from  each  of  7  data  blocks  to  check  the  vari¬ 
ability  of  the  results  one  random  sample  to  another.  Fewer 
than  10%  of  the  averages  differed  by  more  than  1%.  Inner 
products  were  used  to  compare  the  ten  EOFs  that  explained 
the  most  variance  for  a  particular  random  sample  with  those 
based  on  other  random  samples  from  the  same  data  set.  More 
than  half  the  corresponding  inner  products  exceeded  0.98.  A 
value  of  1  indicates  perfect  agreement 

4.  EXAMPLES  OF  EMPIRICALLY  DEFINED 
PATTERNS  OF  SMALL-SCALE 
ATMOSPHERIC  MOTION 

4.1  Averages 

Figures  1  and  2  show  average  patterns  of  small-  and 
medium-scale  variability  for  stable  and  unstable  examples,  re¬ 
spectively.  The  vector  at  the  center  serves  as  a  scale,  pointing 
northeast  at  1  ms“^  The  solid  arrows  indicate  the  horizontal 
wind,  while  the  gray  arrows  show  vertical  motion.  Both 
stable  and  unstable  cases  show  the  average  shear  that  was 
mentioned  earlier.  It  is  stronger  in  the  stable  case.  The  ap¬ 
parent  greater  shear  for  400  meter  resolution  is  because  the 
separation  between  the  planes  is  twice  that  for  the  200  m 
resolution. 


Fig.  1.  Averages  for  two  resolutions  -  stable  case 


Fig.  2.  Averages  for  two  resolutions  -  unstable  case 


4.2  Empiricallv  Defined  Feature 

The  2(X)  m  resolution  EOFs  for  stable  and  unstable  cases 
differ  appreciably,  both  in  their  nature  and  in  the  amount  of 
variance  they  explain.  The  first  two  stable  case  EOFs  (Figure 
3)  explain  37%  and  21%  of  the  variance,  while  the  unstable 
case  EOFs  (Figure  4)  explain  17%  and  12%,  supporting  the 
intuitive  notion  that  an  unstable,  turbulent  atmosphere  dis¬ 
plays  a  much  wider  variety  of  motions  than  a  stable  one. 
Figure  3  shows  tite  dominant  variability  pattern  (EOF  1)  to  be 
a  local  change  in  the  shear.  The  stable  case  EOF  2  also 
changes  the  shear  locally,  but  in  the  opposite  sense  in  the 
lower  and  upper  halves.  Neither  feature  introduces  a  major 
change  in  vertical  motion.  Similar  patterns  were  found  in 
other  stable  cases.  In  one  instance  EOF  1  explained  nearly 
half  the  total  variance. 


Fig.  3. 200  m  resolution  EOF  1  and  EOF  2  -  stable  case 


Fig.  4. 200  m  resolution  EOF  1  and  EOF  2  —  unstable  case 

The  EOFs  for  the  unstable  case  are  characterized  by 
more  vertical  motion  than  for  the  stable.  This  is  not  surpris¬ 
ing  for  a  more  convcctively  active  atmosphere.  They  also 
show  more  evidence  of  rotation.  EOF  1  appears  to  be  a  sec¬ 
tion  of  a  vonex  tilted  about  45®  from  the  vertical,  pointing  in 
the  direction  of  the  shear.  The  vertical  vorticity  component 
would  be  increased  when  the  sense  of  this  feature  is  positive 
and  decreased  for  negative  coefficients.  The  vertical  iix>tion 
is  upward  to  the  northwest  and  downward  toward  the  south¬ 
east  for  a  positive  feature  value.  EOF  2  tends  to  enhance  the 
overall  shear  when  positive.  It  also  introduces  vertical  mo¬ 
tions,  upward  to  the  southwest  and  downward  to  the  north¬ 
east. 
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The  EOFs  determined  frran  the  400  m  resolution 
smoothed  data  are  shown  in  Figures  5  and  6.  Again,  the 
stable  case  EOFs  explain  more  variance  -  40%  and  12%  for 
the  stable  case,  versus  17%  and  14%  for  the  unstable  case. 
EOFl  for  the  stable  case  enhances  the  shear  for  the  400  m 
resoludon  data,  but  mostly  in  die  lower  half  of  the  feature. 
This  reflects  the  nonuniformity  of  shear,  a  tendency  to  de¬ 
crease  with  height  This  EOF  also  has  a  small  component  of 
subsidence  at  several  points,  suggesting  shear  enhancement 
in  subsiding  air.  Values  of  the  400  m  resoludon,  stable  EOF 
2  greater  than  zero  introduce  subsidence  at  the  upper  and 
lower  levels,  enhanced  westerlies  aloft,  and  increased 
southerly  components  below. 

The  400  m  resoludon  EOF  1  for  the  unstable  example 
(Figure  6)  is  similar  to  that  for  the  2(X)  m  resoludon, 
although  the  dlt  of  the  vortex  is  not  as  pronounced.  EOF  2 
for  the  coarser  resoludon  has  litde  effect  on  the  horizontal 
components.  Instead,  it  tends  to  introduce  subsidence  (for 
posidve  values)  throughout  its  volume. 


EOFl  EOF 2 

Fig.  5. 400  m  resolution  EOF  1  and  EOF  2  -  stable  case 


EOFl  EOF 2 

Fig.  6. 400  m  resolution  EOF  1  and  EOF  2  -  unstable  case 


field  with  finer  resolution.  We  expect  to  continue  the  mul- 
tiiesolution  feature  analysis  and  examine  the  fractal  properties 
of  the  flow.  It  win  be  particularly  interesting  to  sec  if  differ¬ 
ent  features  exhibit  different  scaling  properties. 
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